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■ Abstract: We consider the issues that arise out of interpreting the ghost-free bimetric theory 

as a theory of a spin-2 field coupled to gravity. This requires identifying a gravitational metric 
and parameterizing deviations of the resulting theory from general relativity. To this end, we 
first consider the most general bimetric backgrounds for which a massless and a massive spin-2 
fluctuation with Fierz-Pauli mass exist. These backgrounds coincide with solutions in general 
relativity. Based on this, we obtain nonlinear extensions of the massive and massless spin-2 
fields. The background value of the nonlinear massive field parameterizes generic deviations of 
the bimetric theory from GR. It is also shown that the nonlinear massless field does not have 
standard ghost-free matter couplings, and hence cannot represent the gravitational metric. 
However, an appropriate gravitational metric can still be identified in the weak gravity limit. 
Hence in the presence of other neutral spin-2 fields, the weak gravity limit is crucial for 
compatibility with general relativity. We also write down the action in terms of the nonlinear 
massive spin-2 field and obtain its ghost-free couplings to matter. The discussion is then 
generalized to multimetric theories. 
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1 Introduction, motivation and summary 

Theories of interacting spin-2 fields have been considered over many years with various mo- 
tivations, for example, in [1—6], or more recently in [7-15]. Often, these are formulated in 
terms of two metrics, g^ and f^, with non-derivative interactions. These theories generi- 
cally contain Boulware-Deser ghost instabilities [16]. The bimetric theories that avoid this 
problem were written down and proved to be ghost-free in [17, 18]. This was based on [19-21] 
that further developed the massive gravity work in [22, 23], as will be briefly reviewed below. 
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More recently, this was extended to ghost-free theories of many spin-2 fields in terms of J\f 
vielbeins [24], while a formulation in terms of M metrics is given in [25]. 1 

In the ghost-free bimetric theory [17], a priori, the two spin-2 fields g^ and appear 
more or less on the same footing. For obvious reasons, eventually we would like to inter- 
pret this as a theory of a "massive" spin-2 field interacting with gravity. Furthermore, the 
gravity sector of the theory should not show observable deviations from tested aspects of 
general relativity. In this paper we consider such issues that arise when using the bimet- 
ric theory to describe spin-2 field coupled to gravity. The considerations also apply to the 
multivielbein/multimetric case. 

To focus attention, the ghost-free bimetric theory that we will work with has the form, 



S = d 4 x 



m 2 g ^9R(9)+m}^R(f)-2m i ^V(g- 1 f^ n )\ + S m (g,f,^ m ), (1.1) 



with details to be specified later. The particular combination of kinetic and potential terms 
renders the theory ghost-free. The seven parameters of the theory are 77ip, Tfif and five j3 n . 
The simplest possible "matter" interactions that are also known to be ghost-free [17], are of 
the form, 

V=gC g (jg,il;) + yPfC f (f,i/). (1.2) 

Other forms of matter coupling should be explicitly checked for ghost. Generic cosmological 
and localized solutions in this theory could show large deviations from solutions in general 
relativity (GR) although there also exist classes of solutions that are close to GR spacetimes 
[26-32]. 

Below, we will first describe the issues that arise out of interpreting (1.1) as a theory of 
a spin-2 field coupled to gravity, and summarize our results. Then we will briefly review the 
development of spin-2 theories with emphasis on the importance of the nonlinear methods. 

1.1 Issues considered and summary of results 

By construction, in (1.1) around flat backgrounds g^ = = i]^ v , the fluctuations Sg^ and 
Sfuv ar e linear combinations of a massless spin-2 mode 5G^ P (2 polarizations) and a massive 
spin-2 mode 5M^ V (5 polarizations) with a Fierz-Pauli mass term [33, 34]. At nonlinear level 
too, the theory has 7 propagating degrees of freedom, although in that case the analogue of 
the decomposition in terms of mass is not known. 

To regard (1.1) as a theory of a neutral spin-2 field interacting with gravity, one has to 
first identify the gravitational metric, say <7 GR , in terms of g and /. An important restriction 
is that the standard minimal couplings of g GR to matter, demanded by the weak equivalence 
principle, should also be ghost-free. A first guess for g GR , one suggested as far back as [2], 
is the nonlinear extension of the massless mode SG^. But to explicitly check if this allows 



In hindsight, it turns out that Chamseddine, Salam and Strathdee in 1978 [5] had a ghost-free bimetric 
theory, written in terms of vielbeins and with supersymmetry, although the absence of the BD ghost could not 
be demonstrated then. 
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for ghost-free matter couplings, one needs an explicit nonlinear expression for it in terms of 
g and /. The other obvious fall back options are g or /. While not mass eigenstates, these 
have ghost- free matter couplings. 

Having identified a gravitational metric, the next task is to verify that the theory has 
parameter space regions where the solutions for g GR are close enough to GR solutions that the 
bimetric theory is not immediately ruled out on observational grounds. For this, if possible, 
one would like to have some criteria or quantity to parameterize deviations of the bimetric 
theory from GR. In this paper we consider these issues and the results are summarized below. 

Proportional backgrounds and general mass eigenstates: To obtain the mass 
spectrum, we consider the most general class of bimetric backgrounds around which a mas- 
sive mode with a well defined Fierz-Pauli mass term exists. These are the proportional 
backgrounds f^ u = c 2 g^, where c is determined by the parameters of the theory. They coin- 
cide exactly with solutions in general relativity and always exist as bimetric vacuum solutions 
without fixing the parameters of the theory, as long as real solutions for c exist. Flat space 
solutions require fixing only one of the seven parameters by setting the cosmological constant 
to zero. The solutions also exist in the presence of sources, as long as the sources of the g 
and / equations of motion satisfy m 2 T/ lL , = m 2 Tf w . This constraint is not natural, but shows 
that deviations from it drive bimetric solutions away from GR solutions in a generic sense 
(although it is still possible to get isolated GR type solutions). 

Considering fluctuations around fu V = c 2 g^ u backgrounds, one obtains the most general 
expression for the massless mode 5G^ V and the massive mode 5M^ P with a Fierz-Pauli mass. 

Nonlinear modes: We give a procedure to systematically obtain nonlinear combinations 
of / and g that reduce to 5G and 5M at the linear level. Of these we identify one, G^, 
as the nonlinear extension of the massless mode and two possible candidates, M^ v , M^, 
for the nonlinear extension of the massive mode, based on reasonable criteria. These seem 
natural and are simple enough that the expressions relating them to g and / are invertible. 
The vanishing of the nonlinear massive mode, M = 0, is in one-to-one correspondence with 
occurrence of proportional backgrounds = c 2 g^ u . Hence deviations of the VEV of M from 
are driven by the matter couplings of the spin-2 fields and parameterize generic deviations 
of the bimetric theory from GR. If these nonlinear modes have a relevance directly at the 
nonlinear level is not yet answered. 

Identification of gravity: Having a nonlinear massless mode G in hand, we can test 
the conjecture that it should be identified as the gravitational metric. Through an ADM 
analysis we show that within the bimetric framework, the standard minimal coupling of G^ v 
to matter is not ghost-free. This rules out G^ as a candidate for the gravitational metric. 

Another option is g^v (or equivalently f^ u since the formulation is symmetric) which has 
a ghost-free matter coupling. In particular, in the limit m g » cnif, we have 5G — > 5g and 
5M — > 5f. For the nonlinear fields too, in the limit, m g » mf,G^g, although in this case 
M has no particular limit. Hence, if one identifies g^ as the gravitational metric, then in the 
weak gravity limit, g^ v will mostly consist of the massless mode. An obvious consequence is 
that in the presence of massive spin-2 fields, metric perturbations created by a matter source 
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will also have a small massive component. 

Now one can regard the pair g, M as the basic variables and express the bimetric action 
in terms of them. Although the kinetic part in terms of M is more involved than the original 
form in terms of /, the potential is now a finite polynomial in M and does not involve a 
square-root matrix. Also as pointed out earlier, couplings that drive M away from M = 0, 
also drive the solutions for the metric g away from GR. Subsequently, from the ghost-free 
couplings of fav to matter, we obtain couplings between the massive field M and matter. 

Multi spin-2 fields coupled to gravity: Finally we extend the above considerations 
to multi spin-2 theories, as theories of M — 1 massive spin-2 fields coupled to gravity. 

1.2 Background to bimetric theories 

The BD ghost was first observed in massive gravity [16] which corresponds to the bimetric 
theory with one metric held fixed, say, f^ u = r)^. It led to the speculation that such ghost- 
free theories may not exist. The major breakthrough came with the work of de Rham, 
Gabadadze and Tolley [22, 23], who obtained a potentially ghost-free nonlinear massive gravity 
action for = rj^ u , the dRGT model, on which subsequent developments are based. The 
construction was based on a "decoupling limit" analysis (developed for the purpose [12, 
35]) which guaranteed the absence of ghost in that limit. In the perturbative approach, it 
becomes difficult to extend the analysis beyond the decoupling limit, although [23] outlined 
an argument to show the absence of ghost in the Hamiltonian [36] formalism to quartic order 
in = g^v ~ J7/u/- 

To proceed any further, one has to first make sure that the BD ghost is indeed absent 
in the dRGT model for a nonlinear g^ u . One also needs to find out if the more natural case 
of massive gravity with a non-flat reference metric, fn V 7^ T]n V , is ghost-free. It turns out 
that the "decoupling limit", wielded powerfully in [22, 35], is not adequate to address these 
situations, despite some claims to the contrary [37]. 2 

The nonlinear analysis that can answer the above questions was developed in [20, 21] 
based on the formalism of [19]. In [20], for the case — T/fivi h was proved that the dRGT 
model was ghost-free at the nonlinear level. This conclusively established the absence of the 

2 Decoupling limit analysis does not extend beyond decoupling limit for two obvious reasons: (1) The 
decoupling limit involves working with Stiickelberg fields cf> a , introduced via = dfi4> a dv(f> b riab, rather than 
with the metric g^v The <p a mix only with the 4 "gauge" modes of g^v under coordinate transformations 
and learn about the potential BD ghost through them. But the BD ghost mostly resides in the remaining 
6 components of the metric and cannot be completely transferred to the Stiickelberg fields by coordinate 
transformations (otherwise, ghost fluctuations would be expressible as V( M £„) and would not contribute to 
interactions between conserved sources). It is incorrect to argue that the number of independent modes in g^ 
can be reduced to 2 by coordinate transformations in similarity with GR. In GR, this counting is done on-shell 
and holds only for the solutions of the massless Einstein's equations. It is also not possible to transform g^ 
to ?7 M „ by invoking locally flat coordinate systems over the entire spacetime. Obviously, the flat space action 
and equations of motion are not the same as the curved space ones. Thus, away from the decoupling limit 
it is not enough to study the <j> a alone and ignore <7 M „, and the potential ghost within it, by invoking the 
above arguments. (2) So far, it is not obvious how to obtain a decoupling limit for a generic non-flat / M „ . For 
example see [38] for a recent study of the decoupling for a de Sitter / M>/ . 
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BD ghost for the first time. The theory with generic non-flat f^ v was considered first in 
[21] and also proved to be ghost-free nonlinearly. This generic fu V theory provides the most 
natural setup for discussing massive gravity. But from the point of view of this paper, it 
describes rather a massive spin-2 field g^ v in a non-dynamical gravitational background L. 

Finally, [17] obtained the ghost-free bimetric theory for two interacting spin-2 fields g^ w 
and fav with the correct kinetic structure, by exploiting the symmetries of the interactions. 
An issue raised in [39] about the existence of a secondary constraint that was needed for the 
consistency of the formalism was cleared up in [18]. In this paper we will work mostly with 
this bimetric theory. For related work, see [40-47]. 

The rest of the paper is organized as follows. In section 2, we review the ghost-free 
bimetric theory and discuss the proportional background solutions. In section 3, we obtain 
the linear mass eigenstates and mention the weak gravity limit. In section 4, we obtain the 
nonlinear massless and massive spin-2 modes and show that the massless mode does not have 
ghost-free minimal matter couplings. In section 5, we express the bimetric action in terms of 
g and the massive mode M and discuss some of its features. We also discuss the coupling of 
the massive spin-2 field M to matter. In section 6, the discussion is extended to multimetric 
theories. Section 7 contains a brief discussion of the results and some comments. Appendix 
A summarizes some useful equations used in the text. Appendix B describes a rescaling that 
render the action more symmetric. Finally appendix C contains the details of the bimetric 
action in terms of the nonlinear massless and massive modes G — M . 

2 Proportional-background solutions in bimetric theory 

Generic solutions of the bimetric theory have little resemblance to solutions in general rela- 
tivity. In this section we concentrate on a particular class of bimetric background solutions 
that are indistinguishable from backgrounds in general relativity. Although very restrictive, 
this helps in identifying bimetric theories that are close to general relativity. The solutions 
are also useful in analyzing the linear and nonlinear mass spectrum of the bimetric theory. 
We begin with a review of the ghost-free bimetric action. 

2.1 Review of the ghost-free bimetric theory 

Here we briefly review the ghost-free bimetric action and equations of motion. The ghost-free 
bimetric action, excluding matter couplings, is [17], 




) 



(2.1) 



The potential V is given by, 



•l 




(2.2) 



n=0 
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where, e n (X) are elementary symmetric polynomials of the eigenvalues of the matrix X. In 
4-dimensions they can be expressed as, 

e = l, ei = [X], e2 = I([X] 2 -[X 2 ]), e 3 = I([X] 3 - 3[X][X 2 ] + 2[X 3 ]), e 4 = det(X) , (2.3) 

where, [ ] denotes the matrix trace. This potential V was first suggested, for f^ u = rj^ in 
[22, 23] as the unique candidate for a ghost-free massive gravity, based on a "decoupling limit" 
analysis. That it was ghost-free nonlinearly was proven in [20]. The theories with general 
and dynamical where first considered and shown to be ghost-free in [17-19, 21]. The 
square root matrix X = y g" 1 f in V is necessary to avoid the ghost, but also complicates the 
analysis. 

The independent parameters in the action (2.1) are the five dimensionless j3 n and the 
two "Planck masses", m g and mt. The mass scale m is degenerate with the (3 n and can be 
expressed in terms of the other mass parameters. Integrating out matter fields coupled to 
the g and / metrics respectively, results in vacuum energy contributions to fio and /3 4 . The 
remaining f3 n measure the strength of nonlinear interactions between the two metrics. An 
important property of V is, 



V- det g ViVg^f; /3 n ) = V" del / V{y/Th} ; /3 4 -n) • (2.4) 
Then, the action (2.1) is symmetric under the simultaneous replacements, 

g -H / , p n -t /3 4 _„ , m g o m f . (2.5) 
The g^p and J™ equations of motion with generic "matter" couplings are [19], 

RAa) - h,uR(g) - ^ v° u = ^ 7* , (2.6) 

RAf) ~ hUvRU) ~ *r = Tf v . (2.7) 



Here, V$ u and Vj[ u are explicitly given by, 

3 3 

= 9^ Y^iVg^f) , V£, = ^(-l)"/3 4 -n Ux ^(V?^) (2.8) 

n=0 n=0 

where the matrices Y/\ (X) can be expressed as, 

n 

Y "(n)W = E(- 1 ) rXn " re -W- ( 2 - 9 ) 
r=0 

The two expressions in (2.8) are related through the replacements (2.5). In this sense, the 
bimetric theory basically treats g^ and f^ u on the same footing. 

The usual Bianchi identities of the curvature tensors together with jT/^/ = 0, imply 
the Bianchi constraints, which are independent of the scales m g and m/, 

Vf 5) ^ = 0, Vf /) ^ = 0. (2.10) 
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2.2 Proportional background solutions 

Generic solutions of the bimetric theory are very different from solutions in general relativity. 
Here we consider a particular class of bimetric solutions g^ u and f^, sourced by T/L and T^ u , 
which coincide with solutions for the metric in GR. These are solutions of the type 3 

f^ = c% u , (2.11) 

and exist only if t/„ oc Tfa,. This restriction on the matter sources is not always realistic, but 
such solutions are motivated by other considerations discussed at the end of this section. 

For the ansatz (2.11), the Bianchi constraints (2.10) imply that c is a constant. Then 
((2.6) and (2.7)) reduce to two copies of Einstein's equations for the curvatures of g^v, 

- \% V R + A 9 g»u = ±iT$ , iV - \g^R + A/ = , (2.12) 

® f 

where the cosmological constants are given by, 

A 9 = ^ (A) + 3c fa + 3c 2 fa + c 3 /3 3 ) , A/ = ^ (eft. + 3c 2 fa + 3c 3 fa + c 4 /3 4 ) • (2.13) 

& f 

Obviously, the equations are consistent only if, 

(A 9 - A/) g, v = (mf T$ - mf T$) . (2.14) 

The vacuum energy contributions to T^ v and T/L can always be absorbed in fa and fa. 
Hence, the right-hand side can be assumed to contain no piece proportional to g^. Then, for 
localizable sources, each side of the above equation must vanish separately, 4 



' 3 1 



A g = A f , f/, = ^2*. (2.15) 

Later it will be seen that (2.15) is also crucial for the existence of spin-2 massive and massless 
eigenstates. 

The equation A g = Ay determines the constant c in terms of the parameters of the theory, 
m g , uif and the five fa, through the quartic equation, 

a 2 fac A + (3a 2 fa - fa)c 3 + 3(a 2 fa - fa)c 2 + {a 2 fa - 3fa)c - fa = (2.16) 

where 

a = ^l. (2.17) 

m g 

On solving the above equation for c and substituting in A g , one obtains the cosmological 
constant in terms of the parameters of the theory. 



3 For non-proportional metrics, GR type solutions exist for certain choices of j3 n , but for specific metric 
ansatz, say the FRW ansatz [28, 29]. 

4 If A 9 ^ A/, then (2.12) and (2.14) would lead to a complicated differential equation for the T^'s. 
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To get a feeling for the behaviour of c we can solve this equation for the simple case of 
the "minimal" bimetric model, corresponding to (3\ = P% = 0. Then (2.16) gives, 

* = <* = * = ">• (2 ' 18) 

showing that, in general, c 2 could have any value depending on /3o)/?2 and /?4 . This in turn 
gives the cosmological constant, 

= <fc = A = 0). (2.19) 



m 2 3a 2 /?2 — P. 



i 



The scale m 4 can be eliminated in terms of the Fierz-Pauli mass of the massive excitation 
given in the next subsection. 

Note that the most general set of parameters for which the theory admits flat space as 
a solution is obtained from A g = 0, after solving for c. This condition can be solved for one 
of the fi n , leaving the rest free. In contrast, specifying flat space through / = g = r\ and 
A g = Af = will eliminate two of the /3 n leading to a smaller parameter space. For example, 
in the minimal case considered above, theories that admit flat space as a background are 
parameterized by 9/3f = /3o/34, whereas forcing c = 1 gives the smaller parameter space 
A) = Pi = -3/3 3 - 

2.3 Discussion 

If one interprets one of the two metrics, say g^ v , as the gravitational metric coupled to 
ordinary matter T 9 with Planck mass m g , then one recovers all classical backgrounds of GR. 
However, the requirement TjHv = o?T^, v imposed by these solutions is not realistic (except 
possibly for a = 1 such that the two metrics are coupled to the same matter). In spite of this 
such backgrounds are motivated by other considerations. 

(1) The ansatz (2.11) results in the most general class of bimetric backgrounds for which 
there exist a well-defined massive spin-2 fluctuation 5M^ V with a Fierz-Pauli structure, 



m P pV-detff [SMtSM^-(6M^) 2 ] , (2.20) 

along with a decoupled massless spin-2 fluctuation SGuv, as will be discussed below. The 
explicit expressions help extend the linear mass eigenstates 5M and 5G to nonlinear fields M 
and G. The c / 1 case helps in identifying G. Finally, we interpret the bimetric theory as a 
nonlinear theory of a massive spin-2 field coupled to gravity. 

(2) To be consistent with observations, the solutions for the nonlinear field that is iden- 
tified with the gravitational metric must be very close to the corresponding solutions in GR. 
More specifically, corrections to the GR solutions coming from the non-gravitational sector 
must be strongly suppressed in the weak gravity limit. These issues are difficult to investi- 
gate nonlinear ly. Thus as a first step, one can consider perturbations around backgrounds 
of the type (2.11), sourced by independent ST 9 and 5T$ to probe parameter regions that 
suppress deviations from GR. Later we will also see that the nonlinear massive field defined 
with respect to these backgrounds probes deviations from GR. 
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3 Linear massive and massless modes 



In Minkowski backgrounds g = f = r], where the concept of mass is well defined through 
the Poincare group, the spectrum of bimetric theory is known to consist of a massive and a 
massless spin-2 fluctuations [2, 17]. Such backgrounds exist only after two out of the five j3 n 
parameters are fixed. Here we consider the spectrum of linear fluctuations in the theory with 
arbitrary j3 n . In non-flat backgrounds a massive fluctuation is defined as one with a Fierz- 
Pauli mass term (2.20). In bimetric theory, such mass terms arise only around proportional 
backgrounds f^ u = c^g^ considered above. For independent source fluctuations STjh, and 
5T^ V the expressions help in characterizing deviations from GR. The linear mass eigenstates 
are extended to nonlinear fields in the next section. 



3.1 Massive and massless modes in the linearized theory 

Consider canonically normalized fluctuations around the fn V = c 2 g^ u backgrounds, 

9nv = 9nv + ^ &9^v i ffj.v = c g^iy + Sf^ , (3-1) 
Then to linear order, 

( VFW, = c5 p v + 5SP U , where, SS" V = g™ (sf^ - c =g 5g^) . (3.2) 

Expanding V (2.8) and using the results in appendix (A) gives the linearized equations, 

S% Sg P a + A 9 6g^ - ?0 g w (8S" V - 6£6S%) = ±6T$ , (3.3) 

SfZ 5f pa + AfSfc + ^ g w {6S" V - 6£6S%) = ±6T$ , (3.4) 

where, 

S = i(c/3 1 + 2c 2 /3 2 + C 3 /3 3 ). (3.5) 

£ is given in (A. 7). By taking appropriate linear combinations, (3.3) and (3.4) can be easily de- 
coupled in terms of a massive (SM pu ~ g/ixSS x u ) and a massless SG^ ~ Sg^ + cfaf /m g )5f pu 
spin-2 fluctuation. However, this is possible only if A g = Af, which was also required on other 
grounds. Finally, the canonically normalized massless and massive fluctuations become, 5 , 

SGfu, = \ (Sg^ + ca Sf^) , (3.6) 
Vcrcr + 1 

SM »» = i o \ \ ~ ca H^) , (3.7) 

yc z a z + 1 



The canonical normalization is determined from the action requiring that 8g£5g + Sf£Sf = 8ME8M + 
8GSSG. This value will change if £ on the rhs is computed with the background metric G instead of g. 
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where a = — and the notation associates the fluctuations to nonlinear fields G IIV and M uv . 
The corresponding massless and massive equations are, 

E% SG pu + A g 5G, u = bL » v + c ,f V , (3.8) 



m g y/c 2 a 2 + 1 



2 AT (/) - n- 2 /)T (9) 

^ <5M p(T + A fl SM^ + ^ (5M^ - = c ^ f . (3.9) 

mjV cror + 1 

The Fierz-Pauli mass above is parameterized as, 

m| P = m 4 (c/3i + 2c 2 /3 2 + c 3 /3 3 ) f -5— » + A) • ( 3 - 10 ) 

From (3.8) it is evident that, in the background metric g^, the massless fluctuation 5G pu 
couples to matter with the effective Planck mass 



\/c 2 a 2 + 1 = Jm 2 + c 2 m 2 f . (3.11) 



m p = m g 

which must be large for gravity to be weak. This can be achieved in different ways with 
different consequences. It is also evident that at the linear level, SG behaves like the metric 
perturbation in GR. Deviations from GR emerge mainly at the nonlinear level. 

Away from proportional backgrounds, the fluctuations generically do not have a Fierz- 
Pauli mass term. The analysis is further complicated by the fact that in such cases, \/g~ l f 
does not have a simple expansion. 

3.2 Weak gravity limit 

In GR, gravity is described in terms of a massless spin-2 field minimally coupled to matter, 
as required by the weak equivalence principle. Considering the observational evidence in 
support of GR, it is natural to assume that in interacting spin-2 theories too, the gravitational 
interactions must be associated predominantly, if not exclusively, with the massless spin-2 
mode of the theory. The validity of the weak equivalence principle then requires that this 
gravitational mode must couple to matter in more or less the same way that the gravitational 
metric couples to matter in GR. This simple observation leads to the following possibilities. 

(1) Let's assume that the massless mode 5G^ U can be extended to a nonlinear field 
G^y. If Guv could directly couple to matter in a ghost-free manner, using the same minimal 
coupling prescription as in GR, then such matter couplings would not directly violate the weak 
equivalence principle. In this case, one should express the operator £fZ in (3.8) in terms of 
the background G au which will be proportional to g^ u , giving, £fZ{g) = a(c,mf,m g ) £^(G). 
Then the Planck mass is am p which must be large. Later we identify a nonlinear massless 
mode G pu and show that it cannot couple to matter in a ghost-free way. In the absence of 
consistent direct couplings of G au to matter, a different approach is needed. 

(2) Now consider setups where matter fields can directly couple only to the metrics g uu 
or f uu (as in (1.2)), but not to G av . This would be a natural way of accommodating the weak 



- 10 - 



equivalence principle only if or described gravity. On the other hand, empirically, 
gravity is well described by a massless spin-2 field, which in the bimetric setup is G^ u . These 
two requirements can be reconciled if the massless mode is dominated by or fu U . 
Here we consider the possibility that G^ v is mostly made up of g^. 6 The limits in which this 
holds can be identified at the linearized level from (3.6), where, 5G^ U ~ bg^ u holds in the 
limit 



This can be achieved by a small mj or a small c of both. Whether this choice is natural or 
not, will not be addressed here. Also in this limit, the massive fluctuation bM^ is mostly 
saturated by f^. The strength of bM^ interactions depend on the relative values of c and 
nif. Following this reasoning, in section 5 we consider the nonlinear action in terms of g^ u 
and the nonlinear massive mode M^ u . 

4 The nonlinear massless and massive modes 

Now we consider extending the mass eigenstates of linearized bimetric theory to nonlinear 
fields. The ADM analysis of the bimetric action shows that even nonlinearly the theory 
has seven propagating modes [17]. But only their linear fluctuations around = c^g^ 
backgrounds combine into well defined massless and massive spin-2 states. Here we explore 
the nonlinear extensions of these mass eigenstates. In a theory with general covariance, spin- 
2 fields are minimally represented by rank-2 symmetric tensors. Below we find such tensors 
that reduce to the mass eigenstates (3.6) and (3.7) at the linear level. Since this choice is not 
unique, one can also invoke simplicity as a criterion. These are the only criteria employed 
here. We have not considered if the nonlinear modes also propagate two, respectively, five 
degrees of freedom at the nonlinear level. 

4.1 The nonlinear massless spin-2 field G 

The nonlinear massless mode is a symmetric rank-2 tensor G^ p that reproduces the massless 
fluctuation bG (3.6) at the linear level. To determine G^ v we work with the (1, 1) tensor, 



First note that \J g 1 f = g fg l )g- This follows on writing \J g 1 f = y^l + (g 1 f — 1) 
and formally expanding the square-root. We then have the important property, 7 



Now, let us start with a general symmetric (0,2) tensor G^ u (g,f) and, in it, replace / by S 
through / = gS 2 . In general, G^ u could contain powers of 5t, (S" 1 )^, (S T )v and {S~ 1T )v 

6 Equally well, one could replace p M „ by / M „ and c by 1/c 

7 From the above properties of S^ u it follows that, / = S T gS. Hence, S is a local transformation between / M „ 
and g^v, or a generalized vielbein. Further, in terms of S the proportional backgrounds (2.11) are characterized 
by the background value = C(5 M „ which is invariant under general coordinate transformations. 



nig » cmf . 



(3.12) 




(4.1) 



gS = S g. 



(4.2) 
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contracted with g^ u in the right way to produce a (0, 2) tensor. Then using (4.2) it is easy to 
see that general covariance alone restricts G^ u to the form, 

G = g$(S), (4.3) 

where <&t is a matrix function of the matrix St and its inverse, but not its transpose. On 
the proportional background / = c 2 g, where St = c<5t, this becomes, 

G = g$ = g<5>(S = ct) . (4.4) 

Clearly, $ = </>(c)l for a scalar 4>(c). &(S) depends on c in two way: through an explicit 
dependence of $ on c (e.g., through normalizations), and through S. If these two types of 
contributions could be disentangled, $ could be uniquely reconstructed from 5>. 

Let us now consider fluctuations G^ v = G^ u + SG'^ . These can be computed using the 
canonically normalized variables of the previous section. But to ensure explicitly that the 
equations depend on c only through S and not through normalizations, here we work with, 

9^ = 9^ + 6^, 1^ = 1^ + 5/'^, SS^ = ^ x (±5f' Xu -c8g' Xu ), (4.5) 

Then the fluctuation of the nonlinear massless field becomes, 



^ u = 8g'^ x u + g^^- 
c_ d& 



s p 



1 <9$ 



A 



s 



On the other hand, in terms of (4.5) the massless fluctuation (3.6) becomes, 

5G' liy = A{5g' lMU + a 2 5f' llu ) , (4.7) 

with a normalization A. Comparing the coefficients of the fluctuations in (4.6) and (4.7) gives 
two equations for $ as a function of S, evaluated at S = ct, 

A-^ x v = {l + a 2 <?)6\, A' 1 —^ =2ca 2 5 x ji. (4.8) 

OD p el 

The right-hand sides acquire their c-dependence only through S = ct, and not normalizations. 
It is then natural to assume that j4 _1 $(5) depends on c only through S = ct. This leads to 
the unique solution, 

$ = A(t + a 2 S 2 ), (4.9) 
obtained from the first equation on replacing c by S. It gives the nonlinear massless mode, 

G^u = A {g^ u + a 2 f^ u ) . (4.10) 

Without loss of generality we can set A = 1. The fluctuations of this mode can be canonically 
normalized either with respect to g, to give (3.6), or with respect to G^. 
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Other nonlinear extensions of 5G^ U can be found if A 1 <&(S) is allowed to have an explicit 
dependence on c, besides that coming from S. To find these note that, 

6{Sn)Xu =nc n ~ l 5 x Jl. (All) 

cl 



Then general c-dependent solutions of (4.8) can be written as, 

$ = A(a + 2a 2 ^2a n c 2 ~ n S n j , with a + 2a 2 c 2 ^ a n = 1 + a 2 c 2 , ^ n a n = 1 . 



n>0 n>0 n>0 

(4.12) 

Of course, an infinite number of such solutions exists, the one with the lowest power of S being 
= 1 - a 2 c 2 + 2a 2 cS. A non-polynomial solution of (4.8) is A~ l <& = (1 + o^S -2 ) -1 , 
giving the massless mode G" 1 = g^ 1 + a 2 c f . 

Of all these, the c-independent solution (4.10) gives the simplest invertible relation be- 
tween the nonlinear modes and the original bimetric variables g and /. Note that to identify 
this unique c-independent mode, it was important to work with c / 1 backgrounds. Other- 
wise, at c = 1, this criterion is not useful. 

4.2 The nonlinear massive spin-2 field M 

Prom the outset it is evident that the nonlinear massive field is closely related to = 
(-^/ g -1 fY ' ■ This is hinted by the linear equations (3.2), (3.7) and also by the fact that the 
mass potential V (2.2) is a polynomial in S. S^ u is a (1, 1) tensor but it can be brought to 
a symmetric (0, 2) form in more than one way and the nonlinear extensions of the massive 
spin-2 fluctuation are related to these (0, 2) forms. Here we consider two nonlinear extensions, 
M^j, and Mf* u , before discussing the general case. 

In terms of S^ v = g^\S x u , equation (4.2) is the symmetry condition, 

s^ = sv ( 4 - 13 ) 

The fluctuation 5(5^^) = cSg^ + c5M^ v is a mixture of 5g and the massive mode 5M (3.7). 
Now it is obvious that a nonlinear massive mode can be written as, 

M llv = B(g llX S\-cg llv ), (4.14) 

allowing for a normalization B. On proportional backgrounds, M^ u = 0. This is a natural 
vacuum value for a non-gravitational spin-2 field in the sense that it does not break general 
covariance in its vacuum. Fluctuations around this background are the massive modes SM 
(3.7) with a Fierz-Pauli mass term. 8 The condition gS = S T g implies that and g^ u are 
not independent fields, whereas and g^ v can be regarded as independent. 



8 A note on notation: for c = 1 and / M>/ = 77 M „, g~ x M coincides with the K = \/ g^ 1 ^ — 1 in terms of which 
the dRGT model is written. K% was engineered to produce the massive mode around flat space for c = 1, 
while M represents the massive mode around any background for which a Fierz-Pauli mass can be written. 
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A different nonlinear extension of the massive fluctuation is obtained by using G^ v (4.10) 
instead of g^ u , 

M % = A(i+ a ^) ( G ^ sX u ~ cG^) . (4.15) 

The normalization is fixed such that 5M G = 5M. It is easier to invert the relations and 
express (5,/) in terms of (G, M G ) rather than in terms of (G,M). For more on this see 
Appendix C. 

In general, the massive fluctuation 6M has many possible nonlinear extensions. More 
nonlinear extensions can be obtained by following a procedure similar to the massless case. 
By general covariance alone, any matrix function M of g and / can be written as 

M = gV(S). (4.16) 

The fluctuations of this field, 

= v + 9^qs^ 3 SS °P + " " " (417) 
should be equated to the massive fluctuation (3.7) with arbitrary normalization B, 

5M^ v = Bg^5S x v . (4.18) 

This gives, 

S- 1 * = *(d) = 0, B- l —^=8\5^. (4.19) 



dS" p 



s 



Again the c-dependence of the right-hand sides comes only from S and not from normaliza- 
tions. However, now we cannot assume that B~ 1 ^ depends on S and not explicitly on c 
since then ^(c) = would imply ^f(S) = 0, identically. At a least a minimal c-dependence 
is needed to get a nonlinear massive mode with a vanishing background value. The solution 
with the simplest and most natural c dependence is the one corresponding to M x = M (4.14), 



9* V = B [S*-c6*) . (4.20) 

More general solutions of equations (4.19), involving higher powers of S, are given by, 
V = B^2b n c 1 - n S n , with ^6 n , =0, J2nb n = l. (4.21) 

n>0 n>0 n>l 

For example, one can check that the massive mode M G (4.15) is a solution to the above 
equations by reading off the b n from B~ l M G = (1 + c 2 a 2 ) -1 g (— c + S - ca 2 S 2 + a 2 S 3 ). 
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4.3 Absence of ghost-free matter coupling of the massless mode G 

In the previous subsection we obtained a nonlinear generalization G^ v (4.10) of the massless 
fluctuation of bimetric theory. It is natural to ask if could be consistently coupled to 
matter in the standard way, and be interpreted as the gravitational metric. As in GR, such 
matter couplings must be consistent with the weak equivalence principle. Here we show that 
such couplings of G^ u are not ghost-free. Hence regarding G^ u as the gravitational metric is 
not consistent with observations. The alternative, then, is to regard g^ v as the gravitational 
metric and rely on the weak gravity limit discussed earlier. 

Consider standard minimal couplings of G^ u to matter, for example, to a scalar field 0, 

£(m,G) = -V^G G^d^dA. (4.22) 

To see if such ghost-free couplings exist in the bimetric theory, one can perform a Hamiltonian 
(ADM) analysis [36] . We introduce the following notation for the 3 + 1 decomposition of G, 

where Ki = 3 GijKK Standard matter couplings of G^ u such as (4.22), when written in the 
Hamiltonian form using canonically conjugate variables, are linear in K and Ki, 

£ im ,G)=£ + KC + K i C i . (4.24) 

If the dynamics of G^ were described by the Einstein-Hilbert action, \/—GRg ~ n 4 - 7 dt ^G%j + 
KC + KiC 1 , then K and Ki would be Lagrange multipliers in the full theory. Their equa- 
tions of motion would result in four constraints that, along with gauge symmetries, would 
eliminate the ghost and leave two propagating modes for G^. However, in bimetric theory, 
the nonlinear action expressed in terms of G^ v and M (or M G ) is complicated and it is not 
convenient to carry out the ghost analysis in terms of the ADM variables of G^ v . Instead, 
since the bimetric analysis is already known in terms of the ADM variables for g and / [17], 
the strategy here is to analyze the matter coupling (4.24) in terms of these variables. 
To this end, we introduce the 3 + 1 decompositions of g and /, 

The indices on Ni and Li are raised using the inverses of 3 gij and 3 fij, respectively. It is 
known that in terms of new variables n l that parameterize N l — U through (for details and 
the form of the matrix D, see [17, 21]), 

N* = V + Ln* + ND i k n k , (4.26) 

the bimetric theory (2.1) (with no matter couplings) takes the form [17]), 

C = m 2 cS [<K ij d tgij + j'dtfij + NC g + LCf + L0 gf } . (4.27) 
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For convenience we use the scaled fields of appendix B, effectively setting m g = mj. ftij and 
Pij are the momenta conjugate to and fy. The functions C g , Cf and R l g ^ are independent 
of N, L and Lj, but depend on n % and the remaining variables. The action has the additional 
property that the n l equations of motion are independent of N, L and Lj, and determine n l 
in terms of the remaining variables. Thus, N, L and Lj are five Lagrange multipliers whose 
equations, in particular, C g = and Cf = 0, provide the constraints that render the theory 
ghost-free (along with the associated secondary constraints and gauge conditions) [17, 18]. 

To emphasize, this argument for the absence of ghosts crucially depends on the possibility 
of parameterizing the N l in terms of the n l through (4.26). Only then N, L and Lj appear 
linearly in the action (4.27) and enforce the required constraints. Introducing standard matter 
couplings for g^ v and f^ u individually, does not change this story. 9 If, instead, one couples 
a combination of g and / to matter, one has to insure that it does not reintroduce the 
Boulware-Deser ghosts of g and / by destroying the constraints. 

Now, we consider the matter coupling of the nonlinear massless mode G = g + / by 
adding (4.24) to (4.27). The relevant terms in the action are, 

NC g + LCf + Li# gf + +KC + Ki& . (4.28) 

From the 3 + 1 decompositions of G^ u , g^ v and f^ v it is easy to see that, 

Ki = Ni + Li, 3 G 4J = % + % , (4.29) 

K 2 = N 2 + L 2 + s&JKiKj + - V'.V-' - %L l L j . (4.30) 

Already at first glance K is highly nonlinear in N and L which are no longer Lagrange 
multipliers. But this may not yet imply a ghost. Note that after the TV* have been expressed 
in terms of the n l (4.26), we may still carry out a similar reparameterization of the L l in 
terms of some l l . If this could somehow render (4.28) linear in N and L, then the theory 
may still have the constraints to avoid the Boulware-Deser ghosts (although it may propagate 
more than seven modes if the L l constraints are lost). But it turns out that K 2 given above 
is independent of Lj, so reparameterizing it does not help. 

To see this, simplify the expression for K 2 using (4.29) and writing N' 1 = (N l — L l ) + L l , 

K 2 = N 2 + L 2 + (N* - L*){N* - LP) - 3 g ik 3 G kl (4.31) 

Since N l — L l = Lri 1 + ND l kn k , in terms of n*, this is independent of L l and has the form, 

K 2 = ciN 2 + c 2 L 2 + 2c 3 LN . (4.32) 

For the given c±, and C3, this expression is not a prefect square implying that there is 
no way to render K linear in N and L. Hence the associated constraints are lost. Therefore, 
coupling G^y to matter will reintroduce the Boulware-Deser ghosts. Of course, to linear order, 
5 K is linear in SN and SL and, to this order, ghost-free matter couplings exist, as in (3.8). 



9 Coupling g^v and / M „ individually to matter in the standard way, as in (1.2), results in terms of the form 
(4.24), now written for the metrics g and /. Adding these to the bimetric action (4.27) simply modifies C g , Cf 
and Rgf, but keeps the Lagrange multipliers. Hence the no-ghost argument goes through unmodified. 
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4.4 Spin-2 Mixing and Oscillations 

In the absence of ghost-free coupling of the massless mode to matter, one is led to 
consider the standard individual couplings of g^ u and f^ u to matter, which are known to be 
ghost-free. In the weak gravity limit, m g >> nif, we regard g^ u as the gravitational metric. 
The fluctuations 5g and 5f, sourced respectively by ST 9 and ST* , are linear combinations 
of the mass eigenstates 5G and 5M (3.6), (3.7). So the spin-2 states are produced in the 
interaction basis (5g,5f) while they propagate as mass eigenstates (5M,5G). As is well 
known, this will lead to oscillations between (5g, 5f) and a graviton 5g may oscillate to the 
other spin-2 field 5f. This is very similar to neutrino oscillations or the Kq — Kq oscillations. 
The detectability of this effect reduces for higher FP mass of the massive mode. So in cases 
where the massive spin-2 state can be interpreted as a meson or a heavy elementary particle, 
the effect is negligible. But it will have consequences for very light spin-2 states. Nevertheless 
it remains an interesting consequence of the inconsistency of coupling the massless field to 
matter that in the presence of a neutral massive spin-2 field, the gravitational force is mediated 
by a particle that is a superposition of mass eigenstates. 

5 Action for the nonlinear massive spin-2 field 

In this section we consider the bimetric action (2.1) as a theory of a massive spin-2 field M„ v 
in the presence of a gravitational metric g^ u . We also obtain the ghost-free couplings of M^ u 
to fermionic matter fields. 

5.1 The action in terms of g and M 

To regard the bimetric action (2.1) as a theory of a massive spin-2 field in the presence of 
gravity, we express it in terms of the nonlinear massive field and the metric g^. 10 g^v 
couples to matter in the standard way and is the gravitational metric. Mav is a massive 
spin-2 field and couples non-minimally to gravity. The two spin-2 fields mix and their mass 
and interaction eigenstates do not coincide, just as for spin-| fields in the standard model. 
These mixings also result in deviations from GR. In section 3, in the linearized theory, the 
mixings became small in the weak gravity limit m g >> crrif. The hope is that in this limit 
the mixings remain small even nonlinearly and the predictions of this theory do not greatly 
differ from GR. 

To express the bimetric action (2.1) in terms of the fields g and the massive spin-2 field 
M, let us start with the potential V(S ,(3 n ), where S = sf g~ l f (2.2), and use 

S% = g<>°M av + c&%. (5.1) 

10 Other possibilities would be to write the bimetric action in terms of the nonlinear massless field G^ 
(4.10) and the massive field M or M G . The G — M G action is given in the appendix. In the G — M case the 
expressions for g and / become too involved. Since G M „ cannot be coupled to matter in a ghost-free way, here 
we concentrate on the g — M case. 
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Then, the potential written in terms of elementary symmetric polynomials, becomes [19], 

4 4 4 

V(S, /3 n ) = Pne n {S) = Pn c n e n (l + g^M/c) = ^ a c n e n (g^M) = V(g^M, <) 

n=0 n=0 n=0 

(5.2) 

The last step follows from the linear relations between the e n (X) and e n (l + X) . The 
parameters are given in terms of f3 n as 

"4 = & , a c 3 = p 3 + c/3 4 , a c 2 = p 2 + 2c/3 3 + c 2 /3 4 , 

al = Pi + 3c/3 2 + 3c 2 /3 3 + c 3 /3 4 , a§ = /3 + 4c/3i + 6c 2 /3 2 + 4c 3 /3 3 + c 4 /3 4 ■ 

An advantage of writing the theory in terms of M, of course, is that the potential no longer 
involves a square-root matrix. Also, this form is analogous to the familiar form of mass terms 
in field theory, for example, the mass term for massive vector fields, v / — g g fJ,v A fJi A u . 

Now let us turn to the kinetic term for the massive field M which is obtained from 
y/—fR(f) on expressing / in terms of M and g. The criterion is that in the final action for 
M^ u , all covariant derivatives must be with respect to the metric g^ y . To achieve this in a 
systematic way, it is convenient to use 

f/iu = g^p{S Yv i (5-3) 

where S is related to M in a simple way (5.1). Now, the curvatures of / can be expressed in 
terms of curvatures of g using the results in appendix A. In particular, (A. 3) gives, 

iV(/) = RA9) + 2V [M CV - 2C„ { /C aW a , (5.4) 

where V is the covariant derivative compatible with g^ u and 

V = \<f P (V,/^ + Vyffr - Vpfp,) . (5.5) 

Here, / is given by (5.3) and, for ease of notation we have introduced 

<j^ v = {r 1 T V = {S- 2 Y P 9 pV - (5-6) 

Using the curvature relation above, along with y/—f = y/—g det(S) and R(f) = 4>^ u R^ u (f), 
it is a straightforward though tedious exercise to show that (modulo total derivatives), 

v^fR(f) = V^gdet(S) [r u R(g) + n^v^v^J . (5.7) 

The "polarization" tensor II is a function of g and M (through S) given by, 
n^ a/3 = U^^ 5 ^ - 2<t> atl (tP K <t> Sv - ^^W" + ^(pv^As^fy 5^5^ 5 7P S XuJ . (5.8) 

"Explicitly, e (l + X) = e (x), ei(l + X) = 4e (x)+ei(x), e 2 (l + X) = 6e (x)+3ei(x)+e 2 (x), e 3 (l + X) = 
4e (x) + 3ei(x) +2e 2 (x) +e 3 (x), and e 4 (l + x) = e (x) + ei(x) + e 2 (x) +e 3 (x) +e 4 (x). 
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The discarded total derivative terms in (5.7) arise since the left-hand side has fd 2 f terms, 
leading to Md 2 M terms on the right-hand side. These have been converted to dMdM terms 
which amounts to adding the Gibbons-Hawking- York boundary term to the action for /. 

Finally, putting all this together, the action for a massive spin-2 field M interacting with 
a gravitational metric g is given by, 

S gM = f d^xV^g [m 2 g R(g) + m}det(S)^ u R^(g) 

+ m) det(S) n^V^M^V^M^ - 2m 4 V^M, a c n )] . (5.9) 

Here Si, is a function of M^ v . The coupling of M to matter fields will be discussed below. 

The perturbative content of the (g,M) action (2.1) can be discerned easily. The second 
order action for the fluctuations, g = g + 5g and M = + 6M, will contain (Sg) 2 , 5g5M, and 
(5M) 2 terms. It can be diagonalized in terms of the massless mode 5G and massive mode 
5M, leading to the linearized equations (3.8), (3.9). 

5.2 Some features of the g — M action 

Equivalence of the two formulations: It is straightforward that the g and M equations of 
motion obtained from (5.9) imply the bimetric g and / equations of motion and vice versa, 



SSgf | _ SSgM I SM_ _ n SSgf I _ 5Sg M I ^ V I / C 1 

5f l 9~ 5M 'a (5/ ~ ' 5 5 l/~ <J 5 Im (5M ^ 



Hence the g—M formulation is classically equivalent to the g—f formulation and, in particular, 
is also ghost-free. The g — M form has the advantage that it does not involve square-root 
matrices. The price one pays on the other hand, is the tedious kinetic structure for M and 
its kinetic mixing with gravity. Although deriving the equations of motion from (5.9) is not 
convenient, it is much easier to obtain these equations by starting with the g — f equations 
(2. 6), (2. 7), eliminating R(f), and then converting R fiu (f) to R^ u (g) using (5.4). 

In the g—M action, one may perform perturbative calculations even around non-vanishing 
M backgrounds. In the g — f formulation, performing higher order perturbative calculations 
around non-proportional backgrounds is not straightforward as in that case expanding \J g~*f 
is not simple. Of course, the two formulations are not expected to be equivalent in quantum 
theory unless one takes into account the Jacobian factor that arises from the change of 
variables. 

Parameterizing deviations from general relativity: Most of the classical solutions of the 
bimetric action S g f do not coincide with classical solutions in general relativity [26-32] , except 
for the class of proportional backgrounds / = c 2 g considered here, with c determined by the 
parameters of the theory. Generic matter couplings of the g and / metrics will drive the 
solutions away from proportional backgrounds. We are interested in parameter regions where 
these deviations are small. 

In the g — M formulation, proportional backgrounds correspond to M^ u = 0. Hence, in 
the field theory language, a vanishing vacuum expectation value for the massive spin-2 field 
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implies that the classical solutions for the metric coincides with general relativity. Hence, 
deviations of the bimetric theory from general relativity are parameterized by the deviations 
of Muv from zero. These are driven by general couplings of and Mn V to matter fields that 
violate the condition f f = a 2 f 9 (2.15). 

Energy-momentum tensor of spin-2 fields : Consider the gravitational energy momentum 
tensor SS g f /Sguw I n the g — f formulation, fa V contributes to this only through the potential 
V, but its kinetic term y/—fR(f) does not gravitate (of course, it still affects the dynamics of 
g^ v since their equations of motion are coupled). This is while the Noether energy momentum 
tensor computed around flat fuv, will receive contributions from y/—fR(f). 

In the g — M formulation, the gravitational energy momentum tensor SS g M / '^9^v\m 
contains contributions from the kinetic term of M, as well from V and these appear consistent 
with the contribution to the Noether energy momentum tensor around flat g. Of course the 
complete set of equations is the same in both formalisms as these contributions drop out on 
imposing the M equation of motion. The same statements apply to the matter couplings of 
/, that is S(f,ip). In the g — M formulation, M is mostly minimally coupled to g (in the 
sense that the curved space form can be constructed from the flat space expression) except 
for the non-minimal (fi^R^ term. 

Comparison to earlier work : The g — M action (5.9) is useful in comparing to earlier 
attempts of writing a theory of massive spin-2 on a gravitational background. For example, the 
approach in [7, 8] was to start with the quadratic FP theory in flat spacetime and covariantize 
it with a metric g^, also adding non-minimal curvature couplings. This procedure will 
in general not reproduce the action (5.9) as it will miss the factor det(S') as well as the 
complicated polarization structure (5.8), since it only considers terms quadratic in the massive 
field. 

5.3 Coupling massive spin-2 fields to matter 

At present, the only known ghost- free matter couplings in bimetric theory are the standard 
couplings of g^ u and fav to matter sources, as in (1.2). In the weak gravity limit, we inter- 
preted the g^ v couplings as the gravitational interactions of matter fields just as in GR. Then 
the couplings give rise to very specific interactions of the massive spin-2 field with matter, 
dictated by the absence of ghost. To write / in terms of M, again it is convenient to proceed 
through the related matrix S = g^M — cl (5.1) and use, 

ffj,u = S^gafiS^ = g^ a S a ^S^ u , (5-11) 

For fau couplings to bosonic matter, the manipulations are straightforward. For example, for 
a Proca field, £(/, A) = -\\ff f KX F^F vX + 2m\p v A»A V \ , where = 8^A V - d v A^ 
One obtains the coupling of Au to the massive spin-2 field by expressing / in terms of M and 
g through 5. Then, on raising some indices using g^ u one gets, 

C(g,M,A) = - \^ det S [(S~ 2 r p (S~ 2 r a F^ + 2m 2 A (S~ 2 r p A fl AP] . (5.12) 
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The details of this Lagrangian can be investigated further by, e.g., considering flat space 
g = r], and/or expanding S% = c<5^ + ag^ x 8M\ v . 

For fermionic matter a little more work is needed since fermions couple to spin-2 fields 
through vielbeins. For the fn V and vielbeins (below we use the conventions of [48]), 

ffiu = C^t IJabG u , g^u = 6^ i] a i,e v , (5.13) 
equation (5.11) implies the relation, 

e\ = k\e\S\. (5.14) 

Here A is an arbitrary Lorentz transformation, A T r/A = rj. In a Lorentz invariant theory A 
drops out of all expressions, so we set A = 1 without loss of generality. The vielbeins enter the 
couplings through the curved space 7-matrices and through the spin-connections. In terms of 
the Lorentz frame 7-matrices 7 a , one constructs a pair of curved space 7-matrices 7^ = e^a") 0, 
and 7^ = e^al a ■ They satisfy, 

{r,7 b } = 2 v ab , {r,r} = 2/^, {-f,-f} = 2g^. (5.15) 

The inverse of the relation (5.11), with A = 1, then implies, 12 

r = {S- 1 )W'. (5.16) 

Fermions also couple to vielbeins through spin connections in Lorentz covariant deriva- 
tives acting on them, 

D fl = d^-lw; b [j a ,%}. (5.17) 
The spin-connection is given in terms of vielbeins and the Christoffel connection through 

V 6 = + e\ V a % u r^ • (5.18) 

Using (5.14), with A = 1, we can rewrite this as 

= (e^SPMie^S- 1 )-] + e\S\e™{S- l )Zf^) . (5.19) 
From appendix A, V is related to the Christoffel connection V of g pv by 

V = V + , = \r p (V,f pv + V„/ w - V p V) , (5.20) 

where the covariant derivatives are with respect to g^ u , and where f pu is regarded as a 
function of g and M through (5.11). Using these relations, it is straightforward to re-express 
any coupling of f^ u to fermions in terms the massive spin-2 field M^ v and the gravitational 
metric g^ u . The resulting expressions are highly nonlinear in the fields. 

A general A is absorbed by a Lorentz transformation of the spinors, ip' = Aip where, A a b y b = A^j a A. 
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As an example, consider the couplings to lowest order in the fluctuation 5M of the massive 
spin-2 field around the M = background. Then, S^ v = c5u + ag^8M\ v , and to first order, 

V = f <f p [2V [lx 6M u)p - VpSM^] . (5.21) 

Similarly, to this order, 

= &» - %r X 5M x „Y , Dp = Dp - ^[7 p ,7 CT ]V[/M CT]/i . (5.22) 
Let us apply these to the coupling of fp V to a spin- 1 ; field ijj, 

A/2 = i\/^f WTVvl + im^l) + h.c. , (5.23) 

in which ?/5 = ip^j . T>p = + iqAp is the Lorentz and gauge covariant derivative with 
Abelian gauge field A„. We write this to linear order in SM and in the flat space limit 
9pu = Using, y/^f = c 4 + c 3 a5M p p H , one has, 

A/2 =C 3 (l + ^ MP P ) Aree 

- i, (5M^dp + igSMtfAp + mi p ,7 u ]d [p 5M u]tl ) V + h.c. , (5.24) 

where £f ree = i $ {l p dp + iq^Ap + icm^) if). After a partial integration and using (5.15), 
the derivative couplings become, 



—6M U 



L LLV 

cm cS 



Hinder + ^d^-Ud^Y + d> 7 ^ 



+ ^ v d p il)Yi>-\v ilu i>l P dpil) +h.c. (5.25) 

Finally, considering the hermitian conjugate (with the usual hermiticity condition (7^)^ = 
7°7 M 7°), this can be written, 

ic 3 r 1 - - 

5M^ -M^d v +'f& i )\l) -rr*h p d p il) +h.c. . (5.26) 

cm c fr [2 

Couplings of this form were recently considered in a phenomenological context in [49], to 
address the top-quark forward-backward asymmetry. Here, in contrast to [49], the couplings 
are not a priori expected to be flavor violating since they come only from the Lorentz covariant 
derivative and are essentially of a purely gravitational nature. In particular, the first term 
of (5.26) is simply the stress-energy tensor while the second corresponds to a non-derivative 
trace coupling on-shell. 

6 Generalization to more than one massive field 

Recently, in [24] the bimetric action was generalized to a ghost-free theory of N interacting 
spin-2 fields. In this theory, the kinetic term is given in terms of N metrics guv{I) an d the 
interactions between these are constructed in terms of the corresponding vielbeins e a ^(7), 

I d 4 x^gjf) R(I) + ^ I d 4 x U[e(l), • • • , e(M)] . (6.1) 
1=1 J 
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The potential U, constructed in [24], will be presented below in a reformulation. Around flat 
backgrounds, the spectrum consists of one massless and M — 1 massive states. The vielbein 
description is elegant and was very convenient for showing the absence of the Boulware-Deser 
ghosts. For further work, see [25]. 

Here we are interested in interpreting (6.1) as a theory of N — 1 spin-2 fields in the 
presence of gravity. First one has to identify one of the vielbeins, say, e b u (l), with the 
gravitational metric, = e®(l)r] a be b u (l). Then, off shell, the 16 (TV — 1) components of the 
remaining vielbeins contain the 10 (A/" — 1) degrees of freedom for describing N — 1 spin-2 
fields as symmetric rank-2 tensors with kinetic terms consistent with general covariance. In 
addition, there are 6(AA — 1) extra non-dynamical fields, as there are no leftover local Lorentz 
transformations to remove them. The latter, have to be eliminated through their equations 
of motion to isolate the spin-2 content of the theory. A difficulty that arises for N > 2 is in 
disentangling these non-dynamical components from the ones belonging to the spin-2 fields 
in kinetic terms [24]. 

In other words, from the remaining vielbeins, one can construct N — 1 rank-2 tensors 
9^(1) = e^ail) e a v {I) of mixed symmetry. The potential is a function of the #t(7). It is 
difficult to extract from these the spin-2 fields that have kinetic terms, by solving the non- 
dynamical equations. Even more difficult is doing so in a general covariant way. 

This issue is addressed in the metric formulation of the multivielbein action (6.1) that 
was obtained, and argued to remain ghost-free, in [25]. In this setup, the non-dynamical 
fields are isolated from the spin-2 content in a generally covariant way without solving any 
equations of motion, making it appropriate for the considerations here. Then we work with 
the multi spin-2 action, 



J2 [ d'x^gjf) R(l) + ^ / d 4 *r*~* u h . 
i=i J J 



(6.2) 



where Tj 1 i 2 j 3 i 4 , totally symmetric in its indices, contains the free parameters of the theory 
and the multivielbein potential of [24] is reformulated to [25], 



U h ... h = V / -det 5 (l) e^-^e 



x L\{lM g-\l) g{h)] 1 ■■■L\(h)Ug- 1 (l)g(h) 

L v J /Hi L v 



A 4 
M4 



(6.3) 



In this expression the L V X (I) satisfy g^L^p{I)L u a {I) = g pa and L u x (l) 



A' 



They carry 



the 6(AA — 1) non-dynamical parameters as they do not enter the kinetic terms. They can be 
eliminated through their own equations of motion, but solving these equations is not necessary 
to identify the spin-2 content of the theory. In this form, the similarity to bimetric form as 
given in [19] is apparent. 

The equations of motion obtained from (6.1) or (6.2) admit proportional background 
solutions, 



or 



e a M) = c I e! 



(6.4) 
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for I = 2, • • • ,Af, and were we have denoted #^(1) = g^ v and e°„(l) = e a M . The c/ will be 
determined by the parameters of theory. In analogy with the bimetric case, one can introduce, 
for 7 = 2,- •• ,7V, 

S(I) = yjg- x g{l) , M(I) = g S(I) - c ig (6.5) 

Since, 

gS(I) = (S(I)) T , (6.6) 

the Muvil) will be symmetric and represent the N — 1 massive spin-2 fields with vanishing 
expectation values in proportional backgrounds. These generalize the massive mode (4.14) 
of the bimetric theory in the picture that g^ u is the gravitational metric. Note however that 
for generic coefficients Tj 1 i 2 j i j i in (6.2), the actual mass eigenstates will be given by linear 
combinations of the M^ V {I). 

In terms of g and M(I) = gS(I) — cjg, the potential \J—g V[L(I)S(I)) is a finite 
polynomial of its argument. The kinetic terms for the M{I) will simply involve N — 1 copies 
of the corresponding terms in the bimetric case. The L(I) are determined in terms of g and 
the M(J). 

7 Discussion 

The results have already been summarized in section 1 so here we only make some additional 
comments. The nonlinear massless and massive modes were introduced as an extension of 
the corresponding linear modes. It remains to be seen if they have a relevance directly at 
the nonlinear level. Although it is stated that the weak gravity limit is needed to approach 
GR solution in a generic sense, and that the non-vanishing VEV of the massive mode M 
parameterizes deviations from GR, these effects have not yet been quantified. For example, 
note that the Bianchi constraints (2.10) are independent of m g and m,f hence their nontrivial 
consequences will not be affected by the weak gravity limit. 

Another feature of the bimetric theory is that in the g — f formulation with g as the 
gravitational metric, the kinetic energy of / as well as its matter couplings affect gravity only 
through the potential V(g~ l /). In the g — M formulation there are direct couplings between 
g and the kinetic term as well as matter interactions of M. However, after the M equation of 
motion is imposed, the two sectors interact only through V again. In this sense the couplings 
in spin-2 theories are maximally non-minimal. 
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A Curvature relations 



Here we provide a relation between Ricci tensors on a manifold endowed with two covariant 
derivatives (see, for example, [50]). This simplifies manipulation in bimetric theory. 

General relations: Given any two derivative operators V and V, there exist a (1,2) 
tensor field C such that the actions on vectors uo p are related by 

V p u u = V p uj v - C pu a u a . (A.l) 

If V and V are torsion free and compatible with metrics g and g, the tensor C is given by 

C^ a = b aP + V v gpn - Vpg^) . (A.2) 

Defining the associated Riemann tensors by [V M , V a ] uj v = — R pa v^ oj p , it is straightforward 
to derive a relation between the Ricci tensors R pp = R^av* as, 

RA9) = R»M + ^[,c a]u a - ic v{t fc a]p a . (A.3) 

Example: Linearizing General Relativity: Consider a metric g as a perturbation 
around a background metric g, g pv = g~uu + ^9^u ■ To linear order in 5g, the curvature relation 
(A.3) gives 

R^(g) = R,Ag) + 2V [fl 5T a]u a (A.4) 

where, with an obvious change of notation, 5T is given by the linear terms in (A.2), 

5T^ a = \g aP (V M %, + VuSg^ - VpSg^) . (A.5) 

This can be used to expand the Einstein equations, R pv — \g pv R + A.g pu = jjrT pu , to linear 
order in 5g. One gets the background and the fluctuation equations, 



p 



1_ _ _ 1 _ _ l 

Rfiv - -j9nvR + A# Mi , = jpT^ , E p ° v 5g p(T + A5g fM/ = jp ST ^ > ( A - 6 ) 



where we have defined, 



_i 



- g^9 pa V 2 + y>V<V CT - g^BP* + 8ffiR] Sg pa . (A.7) 



Using the background equation, the curvature contributions to £{Z can be re-expressed in 



terms of T pv and A. 



B Bimetric action in scaled variables 



The action (2.1) can be recast in a more symmetric form. Consider the rescalings, 

9 9 / \ TL ~ 9 

m eS ~ t _ m eS t a _( m f\ 5 ™4_ 4™ 



9uu — 9 9uu ) fav — n fay j fin — ( ) fin i m — m o 9 • (B.l 

m z g mj \ m g / m is 
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In the new variables, the action is invariant under the interchanges jf>/,^„H $4-n, 

(B.2) 



S gf = ml s I d 4 * 



gR(g) + J-fR(f)-2m 2 y/-gV[ Vs" 1 /^ 



The analysis of consistency of coupling the nonlinear massless mode to matter in section 4.3 
is performed in terms of these variables. 

C Details of the nonlinear G — M G action 

Below we work with the rescaled variables introduced above. For the sake of completeness, 
here we provide the details for writing the bimetric action in terms of the nonlinear mass- 
less and massive modes G and M G (4.10), (4.15). To do this systematically, note that the 
expressions for G and M G can be inverted to give, 

9„u = G^" 1 )^ = , f„ v = G^- 1 )^ = 4>^ ■ (C.l) 

Here $ a v and are functions of the matrix S given by, 

& v = # + S» a S\ , ^ = SS + {Sr 1 r a (S- 1 )\ , (C.2) 

and S is related to the massive mode M G in a simple way, 

S = G- 1 M G + cl . (C.3) 

These express g and / in terms of G and M G . We also define the inverse matrices, 

4> m Kv = K , ~^ a Kv = SH . (C.4) 

Using (A. 3) the curvatures of g and / can be related to the curvature of G and quantities 
that contain covariant derivatives only with respect to G^ v , 

RAa) = RAG) + 2V { ,c a]u a - 2C u[i fc a] « , (c.s) 

where, 

C^ a = \g afi (V^„ + V v9 fr - Vpgtu,) . (C.6) 

Similarly we have that, 

RAf) = RAG) + 2V^ a - 2C u[ /C a] p a , (C.7) 

where, 

CA = \f aP (V M /^ + V v f fill - Vpfc) . (C.8) 
We further note that the volume densities can be expressed as, 

V 7 - det g = V- det G Vdet S" 1 , and V '- det / = V- det G V det $-! . (C.9) 
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Using these relations we proceed to obtain the general structure of the nonlinear action. 

Kinetic structure for M G : Using (C.l), (C.4), and (C.9), it is a straightforward but 
tedious algebraic exercise to show that (modulo a total derivative), 

V^det^ Q» (2V[,C a] / - 2C v[ /C aV a ) 

= y^fetg" - A U^K\<j> P * + ^ a An P - 25^8^ - r^.ptx?) Va4> Kp V^ Xa ■ (CIO) 

A similar result is obtained for the corresponding term in the f^ u sector by simply replacing 
by $ everywhere in the above (and g by / in the determinant prefactor). Next note that 
from (C.2) we have, 

V a n = (^xS% + S^S\) V a S\ , (C.ll) 

and also 

v a <n = - ((s-'T^s-'Tpis- 1 )^ + (s-'r^s-'r^s- 1 )^) v a s x a . (c.i2) 

These together with (C.10) and its corresponding expression for f^ v give, 

yTdeT^ <T (2V^„ a - 2C v t /C a] /) = y^dit^ V x S\V K S P p , (C.13) 

V^-det/ ^ (2V^„° - 2C„/C tt]/J a ) = V /r deT/ P Q A ^ V A S« CT V^ , (C.14) 

in terms of the polarization tensors and V . We refrain from writing out the full expressions 
for these tensors here and simply note that they can be straightforwardly deduced from (C.10), 
(C.ll) and (C.12). Since 

V K S p p = G? X V K M G , (C.15) 

equations (C.13) and (C.14) provide the kinetic term for massive field M G with only G- 
covariant derivatives. 

Kinetic structure for G: From (C.5) and (C.7) combined with (C.l) we find the 
corresponding relations for the Ricci scalars, 

R(g) = R(G) + S\S a p G pu R^ v {G) + . . . (C.16) 

R(f) = R(G) + (S-y^S-Y^R^G) + ... (C.17) 

where the dots represents the kinetic terms for M G discussed above. Hence, the kinetic 
structure for G is given by the usual Einstein-Hilbert term plus non-minimal coupling of the 
Ricci tensor to S in the g^ sector and to S 1 " 1 in the sector. Apart from this we also take 
into consideration the volume densities given by (C.9). Thus, the full kinetic structure for G 
is given by (omitting the overall factor of y 7 — det G), 

(Vdet S- 1 + Vdet i" 1 ) R(G) + VdetS^ 1 S^S^G^R^G) 

+ 7detl Z i(5- 1 )^(5- 1 ) a /3 G^R^{G) . (C.18) 
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This relation together with (C.13) and (C.14) completely determine the kinetic terms. 

Full nonlinear G—M G action: Using (C.3), the interaction potential is easily expressed 
in terms of M G as in (5.2), 

V 7 - det g V(S, /3 n ) = V-detG (1 + 5 2 )" 1/2 V(M G , a c n ) . (C.19) 

Collecting all the results, we can now write the Lagrangian in (B.2) in terms of G and M G as 

C(G, M G ) = (det(l + S 2 Y 1/2 + det(l + S~ 2 y 1/2 ^j R(G) 

+ (det(l + S 2 )- 1 ^^ + det(l + s- 2 )- l ' 2 V X a l Kp ) V X S\V K S\ 

+ det(l + 5 2 )" 1 / 2 S^S^G^R^G) 

+ det(i + s-'r^is-'r^s^r^R^G) 

- 2m 2 det(l + S 2 )~ 1/2 V{M G , a c n ) , (C.20) 
such that the full action is given by 

S GM = m 2 ff J dW-detG C(G, M G ) . (C.21) 
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